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Abstract 



Acoustic perturbations in a parallel relativistic flow of an inviscid fluid arc con- 
sidered. The general expression for the frequency of the sound waves in a uniformly 
(with zero shear) moving medium is derived. It is shown that relativity evokes a 
difference in the frequencies of the sound-type perturbations propagating along 
and against the current. Besides, it is shown that the perturbations are not purely 
irrotational as they are in nonrelativistic case. For a non-uniformly (with nonzero 
shear) moving fluid a general set of equations, describing the evolution of the acous- 
tic perturbations in relativistic sheared flows, is obtained and analysed when the 
temperature is nonrelativistic. It is shown that, like the nonrelativistic case, in the 
new system: (a) the excitation of vortical, transiently growing perturbations, and 
(b) the excitation of sound-type perturbations, extracting the kinetic energy of the 
background flow, are possible. It is demonstrated that the relativistic character 
of the motion significantly intensifies the efficiency of both these processes. Fi- 
nally, the possible relevance of this study to processes in astrophysical shear flows 
is pointed out. 



1 Introduction 

Acoustic phenomena are abundant in Nature. Compressibility of all known kinds of 
continuous media implies the possibility of excitation and propagation of different kinds 
of acoustic perturbations. In the neutral fluids and gases, we have the plain sound 
waves (Lighthill 1978), in magnetohydrodynamical (MHD) fluids and plasmas, the mag- 
netosonic, and ion-sound waves (Sturrock 1994), and in complex media like the dusty 
plasmas we can excite dust-acoustic and dust ion-acoustic modes (Rao, Shukla & Yu 
1990, Shukla & Silin 1992). 

The nature of all these compressional modes can be seriously affected in differentially 
moving media (shear flows). For example, the perturbations can show large transient 
growths and extract energy from the mean flow (Chagclishvili, Rogava & Segal 1994). 
In a variety of physical situations, where more than one kind of oscillation mode is 
sustainable, the shear flow also allows a mechanism of mutual transformation of the 
modes and hence of mutual exchange of energy between them (Chagelishvili, Rogava & 
Tsiklauri 1996, Rogava, Mahajan & Berczhiani 1996). 

The purpose of this paper is to examine the velocity shear induced acoustic processes 
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for relativistically moving media. This problem, over and above its purely theoretical 
significance, may have applications to a number of astrophysical shear flows [e.g., pulsar 
magnetospheres (Arons & Smith 1979), bipolar and unipolar outflows in Active Galactic 
Nuclei (AGN) (Blandford & Rees 1974)], where relativistic shear flows actually occur. 
We shall "introduce" special- relativistic dynamical effects for the simplest, model parallel 
flow with linear shear; the idea is to single out the principal novelty evoked by the 
relativistic character of the motion. It is hoped that this simple analysis will provide a 
foundation for a more thorough treatment of realistic astrophysical shear flows, which 
are often rather complicated. 

The paper is organized in the following way. In sec. 2 we write down the general 
system of equations of the special-relativistic fluid mechanics. We also consider the 
regular motion of the idealized, plane, parallel "shearing sheet" of the viscous fluid, and 
derive an analytic expression for the mean (equilibrium) velocity profile. 

Section 3 deals with the linear theory of the acoustic perturbations in this flow. We 
first study a uniform (shearless) relativistic flow for which a general expression for the 
frequency of the sound waves is derived. It is shown that when both the mean velocity of 
the flow and the temperature of the fluid are relativistic, there appears the remarkable 
(though expected) difference in the frequencies of the sound waves propagating down 
the current and against the current. Another interesting feature of the sound- type 
perturbations is that they are not curl-free as they are in nonrelativistic hydrodynamics. 
We then introduce shear, and derive the general system of equations describing the 
temporal evolution of the acoustic perturbations. The particular case of the relativistic 
shear flow with nonrelativistic temperature is analyzed in detail. The shear induced 
processes in the flow are rather similar to their nonrelativistic (Chagelishvili, Rogava & 
Segal 1994) analogues: we find almost incompressible, vortex (vortical) perturbations 
that exhibit nonexponential, algebraic transient growth, and the usual, compressible 
(sound-type) perturbations that are able to extract the kinetic energy of the mean flow. 
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It is shown that the rates of the algebraic growth, and also of energy exchange are 
enhanced by a factor 7q (70 being the average Lorentz factor of the flow) resulting in 
tremendous efficiency for extracting energy from the superrelativistic shear flows. 

In the conclusion we briefly discuss the possible areas of astrophysical interest where 
the results of this study can be applied. 

2 Governing equations: mean flow dynamics 

Throughout this paper, we shall use the following notation: (a) The Greek indices will 
range over t, x, y, z and represent space-time coordinates, components, etc.; (b) The 
Latin indices will range over x, y, z and represent coordinates in the three-dimensional 
space. We use geometrical units, so that the speed of light c = 1. 
The spacetime metric is Minkowskian: 

ds 2 = -dr 2 = j] a fidx a dx (i = -dt 2 + dx 2 + dy 2 + dz 2 , (1) 

and the stress-energy tensor for a relativistic, viscous, neutral fluid consisting of particles 
of to , may be written as 

T af3 = (e + P) U a U f3 + Prj *f3 _ 2rj(ya f3 _ £ eh af3 1 ( 2 ) 
in the standard notation. Here U a is the four velocity vector with components: 

LP = [ 7) 7V J ], 7 -(l -VV /2 . (3) 

e is the specific (proper) energy density, P is the specific pressure and 77, £ are the 
coefficients of specific shear and bulk viscosity respectively. The projection tensor h a ^ 
is defined as: 

h af3 =ri a(3 + U a U f >, (4) 
while 6 is the "expansion" of the fluid world lines: 

e=u a , ai (5) 
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and a a @ is the shear tensor: 

<J a p= [U a -^ + Uf,. tlJl h£] /2 - {0/3)h af> . (6) 

Note that a thermodynamic relation between e and P implies: 

e = p+{T-l)- 1 P, (7) 

where p=mon is the proper material density (rest mass density) of the fluid, and n 
is the proper number density of particles. Note also that T — 5/3 for the fluid with 
nonrelativistic temperature and T = 4/3 for the one with ultrarelativistic temperature. 

The equations of motion of the dissipative fluid are contained in the conservation 
law of the stress-energy tensor, 

T4 = 0. (8) 

Let us first consider the equilibrium; we have a relativistic, idealized, plane, parallel 
"shearing sheet" of the viscous fluid. The fluid moves in the x-th direction, its lower 
layer (y = 0) moves with a relativistic velocity Vb<l, while the upper one (y = L) moves 
with Vq + AV (AV<^Vo), L being the width of the "shearing sheet." The symmetry of 
the problem implies that the nonzero component of the fluid mean 3-velocity V x can 
depend only on the y coordinate. The nonzero components of the mean 4-velocity are: 

U* = (1 - ^ 2 )" 1/2 =7, (9a) 
U*=W X , (96) 

Equations (9a) and (9b), coupled with the x-th component of the stress-energy conser- 
vation equation (8) yield: 

d 2 U x , . 

■^- = 0. do) 

with the general solution, U x — Ay + B, which immediately leads to the following 
nonlinear mean velocity profile: 

y/l + (Ay + B) 2 



where A and B are the constants of integration, to be determined through the boundary 
conditions [70 = (1 - V^ 2 )" 1 / 2 ]: 

B = la V , (12a) 
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V + AV 

7o Vo 



yi-(y + A^) 2 

When AF<cVb, the last expression leads to the simple result: 



(126) 



A^I(AV/L). (13) 

^From now on we shall assume that AL<^B, which holds if AVjq<^.Vq- In this case 
(11) leads to the following approximate expressions [a=AV/L]: 

V x (y)~V + ay, (14) 

for the mean flow velocity, and 

7(y)-7o(l + aVb7oV (15) 
for the corresponding Lorentz factor. 

3 Perturbations: linear treatment 

After having established an equilibrium, we now proceed to examine the response of the 
mean flow to linear perturbations. For simplicity, we neglect the viscous terms in the 
stress-energy tensor. The relevant equations are the particle number conservation 

(nU a ), a = 0, (16) 

which in the usual vector notation reads: 

a t (n 7 )+V-(n 7 V) = 0, (17) 

and the equation for momentum conservation. The latter is conveniently derived by 
means of the projection tensor h a/3 , operating on the stress-energy conservation. The 
required equation 

/if T|, a = (e + P)U a U^ a + U,U a P, a + P ti = 0, (18) 



may be written as: 



(e + P) 7 A(7V) + j 2 VD t P + VP = 0, 



(19) 



in the more familiar vector notation. In Eq.(19) D t =[d t + (V, V)] is the convective 
derivative. 

For further analysis it is more convenient to write (17) in the following way 



where g^monj = p"f is the fluid mean density in the laboratory frame. 

Note that due to (15), the mean value of g is spatially inhomogeneous [p = mon]: 



The inhomogencity of the mean flow is the fundamental object of our current inves- 
tigation. Our later calculations show how relativity "creates" enormous effective inho- 
mogeneities as seen in the laboratory frame. We present here a basic physical picture 
of the expected phenomenon: Let us assume that the fluid is inherently homogeneous, 
i.e., its proper mean particle number density n, and its specific mean rest mass density 
p are constants. If we imagine this fluid to be consisting of different "layers," then while 
at rest, the number of "fluid elements" per unit "length" of the layer, as seen in the 
laboratory frame, is the same in all layers. In other words the density n = const. For a 
differentially moving fluid, i.e., when the layers have different velocities in the laboratory 
frame, the situation is much more interesting. Due to the "Lorentz contraction," the 
distances between adjacent fluid particles in each layer (as seen in laboratory frame) 
will decrease; the higher the velocity, the greater the contraction. Thus we arrive at the 
important conclusion: The relativistic differential motion of the fluid induces a density 
stratification (inhomogeneity) in the originally homogeneous fluid. It seems reasonable to 
expect that this special- relativistic effect should play a crucial role in physical processes 
dominated by inhomogeneities. 



D tQ + g (V-V) = 0, 



(20) 



g = p7o(l + aV j$y). 



(21) 
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Let us now decompose the instantaneous values of all physical variables appearing 
in (19) and (20), into their pcrturbational and mean components: 

Q=Q+q', (22a) 

V = V + v, (226) 

7 = 7 + 7= 7o(l + aVo^y + V ^v x ). (22c) 
Using (22a)-(22c), the linearized evolution equations take the form 

[dt + (V Q + ay)d x ]g + p^ {d x v x + d v v v ) + aV^pVy = 0, (23) 



1 

(e+"P)7o 2 ' 

Note that the term aVo^pVy in (23) arises due to the above described relativistic density 



[dt + (V + ay)d x ]v x + av y + - ■ ^ 2 [d t + (V + ay)d x ]P + 4 d x P' = 0, (24) 

7o \ e + -^i7o 

[d t + {V + ay)d x ]v y + 7 —L r - s d y P' = 0. (25) 



stratification effect. 

In writing (23-25), we have neglected the y-dcpcndcnt parts in the coefficients of 
certain terms. However, we do retain the ay term in the convective derivative. This 
term representing shear, even when ay<CVo, is the source of much of the interesting 
physics discussed in this paper. Note that in (23) g is the density perturbation as 
measured in the laboratory frame, while P in (24) and (25) is the specific pressure 
perturbation as measured in the fluid rest frame. It is more convenient to express g 
through the corresponding specific density perturbation in order to have in equations 
all thermodynamic quantities relative to the one and the same frame of reference. The 
result is: 

q = 7p + 7' p~7o(p + pVoilvx). (26) 
Defining d=p jp we can rewrite (23) as 

[d t + (V + ay)d x ] (d + Vo7o w *) + &xV x + d v v y + aVo7o u y = °- ( 27 ) 



8 



The combination (e + P) has a simple physical meaning in either of the extreme 
limits of the nonrelativistic and the ultrarelativistic temperatures. In the former case 
(e + P)—p, while in the latter one (e + P)~4P. In both cases the term is expressible 
through the sound speed C s . For nonrelativistic temperatures C s =(dP/dp) 1 ^ 2 <^l, while 
in the case of the ultrarelativistic temperatures C s = 1/ V3. 

After making the change of variables 



t 1 =t: X\=x- (V + ay)t; y x = y, (28) 

which leads to the following transformation of operators 

[d t + (V + ay)d x ] = d tl ; d x = d Xl ; d y = d Vl - Atxd Xl , (29) 

we can write the closed set of differential equations for three variables: v x , v y and d: 

d tl (d + Vo7o«x) + d Xl v x + (d Vl - Atid Xl )v y + aV Q -ilv y = 0, (30) 

d tl v x + av y + (V Q C 2 s /^)d tl d+ (C 2 J^)d Xl d = 0, (31) 
d tl v y + {C 2 s hl){d m - AtAJd = 0. (32) 
The set of Eqs. (30-32) are now Fourier analyzed in X\ and y\ through the definition 

F = J F(k xl ,k yi ,ti)exp[i(k xl xi + k yi yi)}dk xl dk yi . (33) 

Introducing dimcnsionlcss variables: (r=k Xl ti, R=a/k Xl , (3o=k Vl /k Xl , (3(t) = Pq — Rt), 
we can reduce the Fourier transformed system to the following closed set of first order, 
ordinary differential equations (ODE's) [C 2 =C 2 /^ 2 }: 

d T (d + V 0l %v x ) + i[v x + I3{t)v v ] + RVo^Vy = 0, (34) 

d T v x + Rv y + V C 2 d T d + i(C 2 /^)d = 0, (35) 

d T v y + iC 2 p(r)d = 0. (36) 
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In the nonrclativistic limit (Vb<Cl, 70— 1) these equations reduce to the system that 
was analyzed by Chagelishvili, Rogava & Segal (1994). 

Multiplying (35) by V ^ 2 and combining it with (34), we can write, instead of these 
two equations, a simple pair 



where 5=1 — C 2 V 2 , is limited to the range 2/3 < 5 < 1, and differs noticeably from 
1 only for the case of ultrarelativistic temperatures (C 2 = 1/3) coupled with highly 
relativistic background flow velocity (Vo— 1). In all other cases the parameter is close to 
1. 

From (36) and (38) one can easily deduce 



This relation will be useful in the following considerations. It is also useful to know that 
the parameter R~AV(l x /L)<t.l (£ x is a characteristic length scale of the perturbation 
along X) because both AV and £ X /L are much less than unity. 

3.1 Zero shear (R = 0) case 

We first investigate the shearless (R — 0) system. For this case all coefficients in (36) 
and (38-39) become constants [fi{r) = (3o = const] and the following dispersion relation 



5d T d = i[V C 2 d-(v x +(3v y )], 



(37) 



5d T v x = -SRv y - i(C 2 /rt)d + iVoC 2 (v x + l3v y ). 



(38) 



d T [v y - 6 7 2 pv x ] = 7 2 (SR - iV C 2 p)[v x + pv y }. 



(39) 




sn - V Q C 2 




1 




(40) 



which can be expanded to yield 



n(sn - V0C 2 ) 2 - nc 2 (i/-f 2 + sp 2 ) = o, 



(41) 



gives the normal-mode frequencies, 




(42) 
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Before analysing (42), it is interesting to show that the dispersion relation could be 
obtained by a Lorentz transformation of the sound wave dispersion in the rest frame of 
the fluid, 



^ = d(ki + e u ). (43) 

Now the four-vectors (w, k) in the laboratory frame, and (w, k) in the rest frame are 
the Lorentz transforms of one another, 

w = ^ (w + V k x ), (44a) 

k x =lo{k x - V w), (Mb) 
k y = k y . (44c) 

The sound frequency O, associated with the coordinates (t\, x\, y\), is related to 
the oscillation frequency w by 

w = n-V k x . (45) 

and the corresponding phase is tp=[£lti + k Xl x\ + k Vl yi] = [(Q — Vok x )t + k x x + k y y], 
since k xi — k x and k yi = k y in the shcarless case. 

Inserting (44) and (45) in (43) we get the following quadratic equation 

(1 - ClVt)tf ^-k x Q -%( k -*L + kl)= 0, (46) 
7o 7o \7o / 

whose solution is 

1 



5% 



V C 2 s k x ±C s y/k[+5i$k* ] , (47) 



which is exactly the same as (42); only the normalization is different. 
Let us now examine Eq.(42) in some detail: 

• In the nonrclativistic limit (Vo<Cl, 7o— 1) it reduces to the usual dispersion relation 
describing sound waves: 

fi± =±C s ^/l + Pl (48) 
11 



• for ultrarclativistic temperatures and nonrelativistic mean motion the same rela- 
tion (48) is valid, but now the speed of sound attains its supreme value: C s = 1/ \/3- 

• When the flow mean velocity is rclativistic, but the temperature is nonrelativistic 
the first term in the brackets in (42) is still much smaller than 1 and, therefore, 
the expression for the sound wave frequencies reads: 

n ± = ±(l-V ( ?)C s ^/l + ft/(l-Vi). (49) 

In all these cases the absolute values of the frequencies, describing the sound waves 
propagating along (f2+) and against the background flow, respectively, are equal. 

The degeneracy is, however, destroyed for ultrarelativistic temperatures (when C 2 is not 
too small) and highly relativistic values of the average flow velocity Vq of the shearing 
sheet; the two frequencies, as it is clear from (42), could be vastly different. 

Another interesting feature of the sound waves in a special-relativistic, uniform flow 
is that sound-type perturbations are not purely irrotational (curl-free) as they are in the 
standard nonrelativistic hydrodynamics. We can see this by examining the R — limit 
of Eqs. (34) and (39), 

d T (d + VoioVx) = -i[v x + 0oVy], (50) 
d T [v y - S-fo(3ov x ] = -iVoilC 2 s (3 [v x + f3 v y ], (51) 
which immediately lead to the algebraic relation, 

v y ~ loPo v x = V a C 2 s l3 d + const. (52) 

In the nonrelativistic limit, this relation reduces to the conventional condition v y — 
Pov x — const, which, for sound-type perturbations {const = 0) implies irrotationality of 
these perturbations [ (v y — j3ov x )^{kx v) 2 = 0]. Generally, equation (52), however, docs 
not permit such a conclusion. 
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3.2 Non Zero Shear 

To discuss the nonzero shear case, it is convenient to combine Eqs.(36)-(38) to obtain a 
single differential equation. By straightforward algebra, we can eliminate v x and v y to 
derive an ODE of the third order in the variable d, 

5d*d - 2iV C 2 d 2 T d + C 2 ^ + (3 2 ^j d T d - ARC 2 (3d = 0. (53) 

The other physical variables (v x and v y ) can be determined in terms of d(r), and its 
first and second derivatives: 



n)<)'z<l+ 2V»C 2 O r <l ( ^ +0 2 )<l 



(54) 

v x (t) =i5d T d + V C 2 d- f3v y . (55) 



A full analysis of the problem will require a numerical study of (53), and will be 
attempted at a later stage. To appreciate the new effects introduced by relativity, let us 
consider a strongly relativistic (Vb— 1, 7o^l) shearing sheet with nonrelativistic (Cf<Cl) 
temperature. For this case, the defining equations (37)-(39) can be approximated by 

d T d = -i(v x + pv y ), (56) 

d T v x = -Rv y - i(C 2 h$d. (57) 

and 

d r [v y - l 2 (3v x ] = Rj 2 [v x + (3vy\ , (58) 

It should be noted that the terms proportional to V C 2 /j 2 are negligible even in the 
ultimate case of ultrarelativistic temperatures and superrelativistic mean velocity. This 
coefficient attains its maximum value of 2^/27-0.1283 when V = C s = 1/V3. 

Combining (56)-(58) we can get the following algebraic relation between the physical 
variables: 

v y — "/ 2 /3v x — iR~/ 2 d + const, (59) 
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and derive the following second order inhomogeneous differential equation for v x (t) 

d 2 T v x + %(1 + 7o 2 /? 2 K + constx^ f3 = 0, (60) 
7o 7o 

The general solution of (60) is the obvious sum of its special solution, and the ap- 
propriate solution of the corresponding homogeneous equation: v x = v x + v x . When 
Q(t)=(C s /7q)[1 +7o/? 2 ] 1 ^ 2 nas a s l° w dependence on time r, i.e., 



|r>(r) (1) |«r! 2 (T), 



(61) 



it becomes possible to obtain approximate solution to the homogeneous equation. 

For flows with J?<Cl, condition (61) holds for a wide range of possible values of |/?(t)|, 
implying that the WKB solution 

C 



v x {tY 



exp[i(ip(T) +<p )], 



will be valid for the required range of r. In Eq.(62), the phase 

l p(r) = [ n(r)dr = --§4 \pJl + # + -In 7o /3 + J I + 7o 2 /3 2 



(62) 



(63) 



Due to the smallness of R, the particular solution of (60) can be sought through the 
following series (Magnus, 1976): 



v x (t) = constx ^ R2n Vn{r), 



n=0 



yo(r) = -/3(r)/(l+ 7o 2 /3 2 ), 
Vnir) = -- 



(64) 

(65a) 
(656) 



n 2 (r) dp 2 ' 

Retaining only the lowest order term in R, the approximate solution of the inhomo 
geneous equation (60) may be written as: 



_ C ... , . constx (3{t) 

v x =v x + v x k r __^ exp[i{ip{T) + ifo)\ - 2 



(66) 



' r " /J (l + 7o 2 /? 2 ) 

The value of v x given by (67) can be used to evaluate v y and d. We can then compute 
E{t) — (w 2 + v 2 + d 2 )/2, associated with the specific spatial Fourier harmonic (SFH) 



(67) 
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The energy has two main parts, the first one proportional to Cf2(r) is due to the 
comprcssional perturbations, while the second proportional to const/i}(r) is of a purely 
vortical nature. If the vortical part (mainly incompressible) is dominant, the spectral 
energy of the SFH varies as 

implying transient growth of SFH if f3 — k Vl /k Xl > 0. The energy has a maximum 
for t = t*=(3o/R- This pattern resembles the one obtained in the nonrelativistic case 
(Chagelishvili et al. 1994), but with the remarkable difference that the presence of 7q 
will lead to a much stronger growth for the superrelativistic case. In particular, the ratio 
of the maximum value of the SFH spectral energy (attained at the moment r = ) to 
its initial value is: 

l+7o/3o- (69) 



E m in E(0) 

which, unlike its nonrelativistic counterpart, shows substantial increase even if f3o<l. 
Thus relativity in conjunction with shear can induce large transient growths for a wide 
range of k Vl and k Xl . Significant growth will be seen in regimes where there was virtually 
no growth in the nonrelativistic case. And in the regimes of nonrelativistic growth, the 
relativistic growth is Jq times larger. 

For the compressional perturbations, (C^O, const = 0) the energy change takes the 
form 

S(r)~[l + 7o 2 (A)-i2r) 2 ] 1/2 . (70) 

If k yi and k Xl are such that (3q > 0, then for < t < r*, the energy decreases and 
reaches its minimum at r = r*. For r* < t < oo, the SFH "emerges" into the area of 
k-space in which k y (T)k Xl < 0, and begins to grow. If initially (3$ < 0, the perturbation 
energy monotonically increases. The perturbations grow up effectively extracting the 
energy from the kinetic energy of the mean (background) flow. 

Like the case with vortical perturbations, the presence of Jq in (70) intensifies the 
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process of the extraction of the mean flow energy by the sound-type perturbations. The 
rate of energy increase in the "growth area" (characterized by the derivative d T E) is 
proportional to 7q, and is quite fast for large values of the Lorentz factor characterizing 
the mean flow. 

4 Conclusion 

In this paper we have investigated the novel features induced by relativity in the acoustic 
phenomena by examining a simple, model parallel flow with linear shear. The acoustic 
perturbations are excited in a background flow corresponding to an idealized, plane, 
"shearing sheet" of a viscous fluid moving with relativistic speed. Principal results of 
this study are: 

In the shearless limit, when both the temperature of the medium and the mean 
velocity of the flow are strongly relativistic, we find a noteworthy difference in the 
frequencies of the sound waves propagating along and against the mean flow. Relativity 
also causes the violation of the irrotationality of sound- type perturbations. 

For the sheared system, we have first derived the general system of equations for 
the temporal evolution of the acoustic perturbations. We then analysed in detail the 
particular case of the relativistic shear flow with a nonrelativistic temperature. It was 
shown that the standard shear-induced processes like the nonexponential, algebraic tem- 
poral evolution with transient growth of the vortical perturbations, and the ability of the 
compressible (sound-type) perturbations to extract energy from the mean flow, persist 
in relativistic dynamics. In fact, these processes become much more efficient; the rate of 
a typical process is enhanced by the square of the Lorentz factor 7g which can be very 
large for superrelativistic flows. 

It is likely that the model flows discussed in this paper will have common features 
with real relativistic astrophysical flows. The possible candidates are the flows in the 
relativistic regions of accretion discs in quasars, jet-like outflows associated with Ac- 
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tive Galactic Nuclei (Granik & Chapline 1996), pulsar winds, and flows in the inner 
magnetospheres of pulsars (Lyubarsky 1995). In many of these real astrophysical flows 
the special-relativistic effects will be intricately interlaced with the general-relativistic 
effects evoked by strong gravitational fields, and also with the electromagnetic effects 
(for plasma flows) arising due to the existence of superstrong magnetic fields. We hope 
that the methods used and the results obtained in the present study will serve as start- 
ing points for some of the future investigations of these interesting and complicated 
relativistic flows. 
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